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^ Abstract 

Leptogenesis with heavy neutrino flavours is discussed within a density matrix 
i_i formahsm. We write the density matrix equation that describes the generation of 

^ the matter-antimatter asymmetry, for an arbitrary choice of the right-handed (RH) 

^ neutrino masses. For hierarchical RH neutrino masses lying in the fully flavoured 

regimes, the density matrix equation reduces to multiple-stage Boltzmann equa- 
^1 tions. In this case we recover and extend results previously derived within a quan- 

ts;! tum state collapse description. We confirm the generic existence of phantom terms, 

which are not washed out at production and contribute to the flavoured asymmetries 
proportionally to the initial RH neutrino abundances. Even in the A'^i -dominated 
^ scenario they can give rise to lepton flavour asymmetries much larger than the 

^ baryon asymmetry with potential applications. We also conflrm that there is a 

^ (orthogonal) component in the asymmetry produced by the heavier RH neutrinos 

which completely escapes the washout from the lighter RH neutrinos and show that 
phantom terms additionally contribute to it. The other (parallel) component is 
washed out with the usual exponential factor, even for a weak washout. Finally, 
as an illustration, we study the two RH neutrino model in the light of the above 
findings, showing that phantom terms can contribute in this case as well. 



1 Introduction 



Leptogenesis |lj is a direct cosmological application of the see-saw mechanism [2j for the 
explanation of neutrino masses and mixing and it realises a highly non trivial link between 
cosmology and neutrino physics. The discovery of neutrino masses and mixing in neutrino 
oscillation experiments [3] has drawn great attention on leptogenesis that became the most 
attractive model of baryogenesis for the explanation of the matter-antimatter asymmetry 
of the Universe. 

In most cases, classical Boltzmann equations are sufficient for the calculation of the 
final asymmetry |H El El El IS]- However, when lepton flavour effects are taken into 
account [IHl HI] , different set of classical Boltzmann equations apply depending whether 
the asymmetry is generated in the one-flavour regime, when the mass of the decaying RH 
neutrinos Mj is much above 10^^ GeV, in the two-flavour regime, for 10^^ GeV ^ Mi ^ 
10^ GeV, or in the three-flavour regime for Mj <^ 10^ GeV. Moreover classical Boltzmann 
equations fail in reproducing the correct result in the transition regimes for Mj ~ 10^ GeV 
and for Mj ~ 10^^ GeV. However, in the case that just the lightest RH neutrino species 
is assumed to be responsible for the generation of the asymmetry, classical Boltzmann 
equations provide quite a convenient description in phenomenological investigations, since 
the flnal asymmetry can be expressed in terms of simple analytical expressions that well 
approximate the numerical solutions [9l [T2] . 

On the other hand, the contribution from heavier RH neutrinos can also relevantly 
contribute to the flnal asymmetry (heavy neutrino flavour effects) and has, therefore, 
consistently to be taken into account in general [13]. When lepton flavour effects are 
considered [2], a reliable calculation of the asymmetry cannot neglect the contribution 
from the heavier RH neutrinos even in the two RH neutrino model [15] usually considered 
as a paradigmatic case for the validity of the traditional A'^i-dominated scenario. It has also 
been shown that a successful A^2-dominated scenario is naturally realised in an interesting 
class of 5*0(10) inspired models [TB] . 

When heavier RH neutrinos are included, one has to distinguish quite a large number of 
possible mass patterns with different corresponding sets of classical Boltzmann equations 
for the calculation of the flnal asymmetry. For example, in the typical case of three RH 
neutrinos one has ten different possible mass patterns [IT] (see Fig. 1). In addition, the 
requirement that all RH neutrino masses do not fall in a transition regime becomes clearly 
much more restrictive. 

Moreover new effects arise when heavy neutrino flavours are taken into account. First, 
part of the asymmetry generated by a heavier RH neutrino, the flavour orthogonal com- 
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Figure 1: The ten different three RH neutrino mass patterns requiring 10 different sets of 
Boltzmann equations for the calculation of the asymmetry jl7] . 



ponent, escapes the washout from a lighter RH neutrino |B]. Second, parts of the flavour 
asymmetries (phantom terms) produced in the one or two flavour regimes are not washed 
out at the production [18]. 

Therefore, it is necessary to extend the density matrix formalism beyond the tradi- 
tional A^i-dominated scenario pT| [TU] and account for heavy neutrino flavours effects in 
order to calculate the final asymmetry for an arbitrary choice of the RH neutrino masses. 
This is the main objective of this paper. At the same time we want to show how Boltz- 
mann equations can be recovered from the density matrix equations for the hierarchical 
RH neutrino mass patterns shown in Fig. 1 allowing an explicit analytic calculation of 
the final asymmetry. In this way we will confirm and extend results that were obtained 
within a simple quantum state collapse description. For illustrative purposes, we proceed 
in a modular way, first discussing the specific effects in isolation within simplified cases 
and then discussing the most general case that include all effects. The paper is organised 
in the following way. 

In Section 2 we discuss the derivation of the kinetic equations for the A^i-dominated 
scenario in the absence of heavy neutrino flavours. This is useful both to show the 
extension from classical Boltzmann to density matrix equations and to highlight some 
features that will prove to be quite important when, in the next Section, we include 
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heavy neutrino flavour effects. In particular we show the existence of phantom terms and 
how the expression for the CP asymmetry matrix can be unambiguously derived from the 
flavoured CP asymmetries, taking into account the different flavour compositions of lepton 
and anti-lepton quantum states. We also interestingly point out that, though phantom 
terms cannot affect the final asymmetry in the A^i-dominated scenario, they can give rise 
to lepton asymmetries much larger than the observed baryon asymmetry with potential 
applications. In Section 3 we start discussing the case where two heavy RH neutrino 
flavours are involved directly in the generation of the asymmetry, considering a simplified 
case with only two charged lepton flavours. In this way we simplify the notation and we 
better highlight the main results. In this Section we are particularly interested to show 
two effects that specifically arise when the interplay between heavy neutrino and lepton 
neutrino flavours is considered. The first one is phantom leptogenesis [18]. The second, 
that we call projection effect, is how part of the asymmetry generated by a heavy RH 
neutrino, the component orthogonal to the heavy neutrino flavour associated to a lighter 
RH neutrino, is not washed out by the inverse processes of the latter 120] ■ We also show 
that these two effects, phantom leptogenesis and projection effect, in general combine with 
each other. In Section 4 we extend the discussion to the general case with three heavy 
neutrino flavours and three charged lepton flavours. In this section we finally obtain 
general density matrix equations for the calculation of the asymmetry for an arbitrary 
choice of the RH neutrino masses. From these equations we derive the classical Boltzmann 
equations for a particularly interesting case: the two RH neutrino model. The derivation 
can be easily extended to all 10 hierarchical RH neutrino mass patterns shown in Fig. 1. 
In Section 5 we draw the conclusions. 



2 Kinetic equations for the A^i dominated scenario 

We discuss leptogenesis within a minimal type I seesaw mechanism with three RH neutrino 
species, Ni, N2 and N3, with masses Mi < M2 < M3 respectively. 

In this section we review the main steps underlying the derivation of the kinetic equa- 
tions in leptogenesis when heavy neutrino flavours are neglected, assuming that only the 
lightest RH neutrino decays and inverse processes contribute to the final asymmetry: the 
traditional A^i-dominated scenario. 

We first derive the Boltzmann (rate) equations and then we extend them writing the 
density matrix equations, accounting for quantum decoherence and flavour oscillations. 
This derivation will prove to be useful not only to setup the notation but also to highlight 
some basic features of the kinetic equations in leptogenesis that will be relevant when we 

4 



will include heavy neutrino flavour effects in the next section. 

We will neglect different effects, processes and corrections that have been studied 
during the last years and that will not play a relevant role in our discussion. These include 
for example AL = 2 washout [H E], AL = 1 scatterings [5], momentum dependence [2l] . 
thermal corrections |8l |22] , flavour coupling from the Higgs and quark asymmetries |6l [18] , 
quantum kinetic effects [23] • 

We will moreover always assume vanishing initial asymmetry though notice that the 
results that we will obtain in Section 3 when we discuss the projection effect, are also 
important in order to describe the evolution of a non-vanishing pre-existing asymmetry 
[201 [T7|. 



2.1 Boltzmann equations 

If we indicate with Fi the decay rate of the lightest RH neutrinos into leptons, A^^i — )■ 
ii + and with Fi the decay rate into anti-leptons, A^i — )■ £i + $, we can introduce the 
decay term Di and the washout term Wi given respectively by 

F -I- F / 1 \ 1 -I- 1 

D^{z) = \~^^ =KiZ (-) and Wi{z) = - \~[ ^ = - /Ci(z) , (1) 



Hz ^ ' 2 Hz 4 

where z = Mi/T, Ki = (Fi + Ti)t=o/Ht=Mi is the decay parameter, H is the expan- 
sion rate and the averaged dilution factor, in terms of the Bessel functions, is given by 
(1/71) = )Ci{z) / )C2{z) . Under the afore mentioned assumptions and approximations and 
considering the unflavoured regime, the calculation of the asymmetry is described by the 
most traditional set of kinetic equations for leptogenesis from the decays of the lightest 
RH neutrinos A^i P [9] 



dz 
dNs-L 



-Di (iV^, - iV^^J , (2) 

eiD,{NN,-N'^\)-W,NB-L, (3) 

where with Nx we indicate any particle number or asymmetry X calculated in a portion 
of co-moving volume containing one heavy neutrino in ultra-relativistic thermal equilib- 
rium in a way that N'^_[T ^ Mi) = 1. In this way baryon-to-photon number ratio at 
recombination is related to the final B — L asymmetry by 



r]B = asph^^0.01iV],_^, (4) 

7 

to be compared with the value measured from the CMB anisotropics observations [2^ 

ryg^'^ = (6.2±0.15) X 10-^°. (5) 
5 



Let us very shortly recall the basic steps for the derivation of the Eq. ^ for the B — L 
asymmetry. Ignoring the reprocessing action of sphalerons, we can write 



(6) 



dz dz dz 

The net production rate of leptons and anti-leptons is then given by the difference between 
the production rate due to decays and the depletion rate due to inverse decays, for leptons 

^^^^ N ^'"^ N (7) 

and for anti-leptons 

-ir-w-z^'^^'irz^'^- 

The inverse decay rates are related to the decay rates by [5j 

r(^ = ri^ and rf = r,^, (9) 

where N^^ = N^'^ = 1 is the number of leptons ii and of anti-leptons ii in thermal 
equilibrium for vanishing asymmetry. The number of leptons and anti-leptons can be 
then recast as 

Ni, = I {Ne, + iV,-J + ^ (iV,, - iV,-J (10) 



and 



Ni, = I {Ne, + iV,-J - i {Ne, - iV,-J (11) 

Inserting these last expressions into the Eq. ([6| and neglecting terms OlN^^^i^), the Eq. (|3 
is obtained, with Di and Wi given by the Eqs. 

The solution for the final asymmetry has a very simple analytical expression [9J 



Nb_i = e\ k{Ki) , with k{x) 



X Zb[X) 



1 - exp ( --xzb{x^ 



(12) 



^This expression directly accounts for the resonant AL = 2 contribution that is needed not to violate 
the Sakharov condition on the necessity of a departure from thermal equilibrium for the generation of 
an asymmetry jjj. Here we are just interested in showing the separate Boltzmann equations for lepton 
anti-lepton numbers that we will use in the next subsection to derive the CP violating term in the density 
matrix equation. 
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where k{Ki) is the final efficiency factor that here we have written in the simple case of 
initial thermal A'^i-abundance. The asymmetry is generated in a quite narrow interval of 
temperatures around Tbi = Mi/zsi, where zbi = zb{Ki) [9J. 

The unflavoured assumption, underlying the Eqs. ^ and ([s]), proves to describe the 
correct final asymmetry only for masses Mi > 10^^ GeV [25]. In this range of masses, 
during all the relevant period of the asymmetry production, the lepton and anti-lepton 
quantum states produced from the decays of the A^i, that we will indicate respectively 
simply with |1) and can be treated, in fiavour space, as pure states between their 
production at decay and their absorption at a subsequent inverse decay. They can be 
expressed as a linear combination of fiavour eigenstates (a = e, fi, r) 

|l) = ^Ci„|a), Cia = {a\l) and |1) = ^ Ci„ |a) , Ci^ = , (13) 

a a 

where notice that in general the final anti-lepton states produced by the N2 decays are 
not the CP conjugated of the final lepton states and therefore, in general, Cia 7^ Cia [TO] . 

It will prove useful to introduce the branching ratios pia = jCiap and pia = 
giving respectively the probabilities that a lepton ii or an anti-lepton ii is found in a 
fiavour eigenstate a or a in a fiavour measurement process. It is also useful to recast the 
branching ratios as 

Pi«=< + ^, Pi«=<-^, (14) 

where the average p^^ = (pia+Pia)/'^ are the tree level contributions, while the differences 
^Pia = Pia — Pia describe the CP violating contributions from the interference with loop 
diagrams. 

If the unffavoured regime holds during the period of generation of the asymmetry, for 
z < zb, the ffavour compositions of the leptons and of the anti leptons do not play any role 
since the only relevant interactions, the neutrino Yukawa interactions, are ffavour blind. 
On the other hand for masses 10^^ GeV ^ Mi ^ 10^ GeV, the coherent evolution of the 
|1) and |1) quantum states breaks down before they inverse decay interacting with the 
Higgs bosons, due to collisions with right-handed tauons. At the inverse decays, lepton 
quantum states can then be described as an incoherent mixture of tauon eigenstates |r) 
and of \t^) quantum states. These second ones are a coherent superposition of muon and 
electron eigenstates that can be regarded as the projection of the lepton quantum states 
|1) on the plane orthogonal to the tauon ffavour (see Fig. 2). In this two fully ffavored 
regime, classical Boltzmann equations can be still used as in the unffavored regime, with 
the difference, in general, that now the ffavour compositions of leptons and anti-leptons do 
play a role in the generation of the asymmetry. In this case each single ffavour asymmetry 

7 



Figure 2: For 10^^ GeV ^ Mi 3> lO^GeV, the lepton quantum states |1) can be treated 
as an incoherent mixture of a r and of a component during the generation of the 
asymmetry and a two fully flavoured regime applies. 



has to be tracked independently and the total final B — L asymmetry has to be calculated 
after freeze-out as the sum of the two flavoured asymmetries, a r asymmetry and a 
asymmetry. To this extent, we have to introduce the flavoured CP asymmetries 



Pia 



Pia 



Pia 



^P^. 



(15) 



Since sphaleron processes conserve the quantities = B/3 — {a = e, /i,r) [H], these 
are the convenient independent variables to be used in the set of Boltzmann equations 
that can be written as 



dz 
dN^^ 

dz 
dNA . 



Sir Di {Nn, - iV^^J - p?, Wi Na^ 



(16) 



dz 



where p'j'^x = Pie + P% and e^^x = eie + £1^. 

Using the decomposition of the flavoured CP asymmetries in terms of p^^, and Apia 
(cf. Eq. (15)) and assuming strong washout for both flavours (i.e. Kir^Ki^i^ ^ 1), the 
final asymmetry is approximated by the expression 



2eiK{Ki) + 



(17) 



where Kia = p^^ ^i- This approximated expression shows how, compared to the expres- 
sion obtained in the unflavoured case, large lepton flavour effects can arise only when 
leptons and anti-leptons have a different flavour composition, for non-vanishing Apir ^ 
The most extreme case is realized when ei = ^U\. In the unflavoured case this would 
imply a vanishing final asymmetry but in the flavoured case it does not. It should be 
indeed noticed that when flavour effects are considered, B — L violation is not a necessary 
condition for the generation of a baryon asymmetry via leptogenesis, it is sufficient to have 
a A„ violation accompanied by an asymmetric washout between the two flavours, that in 
this context corresponds to the requirement of departure from thermal equilibrium. 

For Ml <^ 10^ GeV muon interactions are able to break the coherent evolution also 
of the \t^) quantum states between decays and inverse decays during the period of the 
generation of the asymmetry. In this case a three-fuUy flavoured regime is realised and 
the set of classical Boltzmann equations is a straightforward generalisation of that one 
written in the two fully flavoured regime. In the iVi-dominated scenario with hierarchical 
RH neutrinos, because of the lower bound Mi > lO^GeV for successful leptogenesis [2S1E], 
a three fully flavoured regime is not relevant for the calculation of the final asymmetry. On 
the other hand, in a iV2-dominated scenario, a three flavoured regime has to be considered 
in the calculation of the lightest RH neutrino washout [H] . 

2.2 Density matrix equations 

Within a density matrix formalism [HI [271 HH], the description of leptogenesis is more 
general than with classical Boltzmann equations, since it makes possible to describe those 
intermediate regimes where lepton quantum states interact with the thermal bath via 
charged lepton interactions between decays and inverse decays but not so efficiently that 
a quantum collapse approximation can be applied in a statistical description. In this case 
the lepton quantum states cannot be described neither in terms of pure states nor as an 
incoherent mixture. Yukawa interactions and charged lepton interactions compete with 
each other in the determination of the lepton flavour. A statistical quantum-mechanical 
description of lepton flavour cannot treat leptons as decoupled from the thermal bath. 
Therefore, the concept of lepton quantum states itself is blurred since one should con- 
sistently describe together leptons and thermal bath. A density matrix formalism [28] is 
then particularly convenient since it still allows to describe the leptonic subsystem in a 

^Notice that relaxing the assumption of strong washout for both flavours one can only get an asymme- 
try that is even closer to the unflavoured calculation. Indeed in the limit of no washout {Kit, ^irj^ 1) 
one exactly recovers the unflavoured expression for the final asymmetry. 
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separate way, neglecting back-reaction effects and encoding tlie coupling with the thermal 
bath in the evolution of the off-diagonal terms of the lepton density matrices. 

Let us see how a density matrix equation for the B — L asymmetry can be obtained 
starting first from the case where charged lepton interactions are negligible. In this case 
we just expect to reproduce the Eq. ([3]). 

Let us consider a simple two lepton flavour case able to describe the intermediate 
regime between the unfiavoured case and the two fully flavoured regime that are recovered 
as asymptotic limits. The two relevant flavours are then r and (see Fig. 2). In this 
two flavour space the flavour composition of the lepton quantum states produced by the 
A''i decays can be written as (a = r, rf*-) 

|1) = Ci.|r)+Ci,x|ri^), Ci„ = («|l), (18) 
|i) = Ci,|f)+Ci,x|fi^), Ci^ = (a|l). (19) 

At three level, the amplitudes Cia and Cja are given by 

Cl= /"^ and . (20) 

Including one- loop CP- violating corrections, these amplitudes become 

^ v/(/^t/,),- 2 Re(/.t + ^''^ 

C^a = ^ {hai + {h''^v)ai) . (22) 



We are following here notation and formalism introduced in [29] and more recently in [30] . 
The one-loop corrections are included in the ^„ and ^„ functions, which are given by 

[UM^)],^ ^ [u^{Mf) + Mb{M^){h^hrM]^^, (23) 
[UM^)],^ ^ [v'^{M^) + Mb{M^){h^h)M]^^. (24) 

Note that the mass matrix being diagonal, we simply have Mki = MiSki- The flrst term 
on the right-hand side describes the self-energy correction, whereas the second one is the 
vertex correction. The imaginary part of the loop factor b{M^) evaluated on mass shell 
for the RH neutrino Ni is given by 



where f{x) = y/x (l — (1 + x) log {■^^))- As for the u and v terms in Eq. (23), they are 
given by 

u^k{M^) = u.kiM^) [MiE^MM-) + MkJ^NMiM^)] , (26) 
v,k{M^) = oo,k{M^) [M,EN,HiM!) + Mk^NAM-)] . (27) 
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They depend on the self-energy T,iy^ki{Mf) = a(Mf){h^ h)ki, where the loop factor a, 
whose imaginary part is Im(a) = — l/(167r), is evaluated again on mass-shell for the RH 
neutrino Ni. Finally, the propagator u, evaluated on shell, can be found to be 



In the two flavour bases li-li and li-^i-, the lepton and anti-lepton density matrix 
are respectively simply given by plj = diag(l,0) and p|j = diag(l,0), where i,j = 1, l"*" 
and i,j = 1, l"*". Because of the different flavour composition of leptons and anti-leptons, 
the two bases are not CP conjugated of each other 

If we introduce the lepton and anti-lepton number density matrices, respectively N^^ = 
Ne-^ and N-j = N^^ pi-., their evolution aX T Tl is given by 

^^^1 - TFT ^^^1 PiJ ' -IT - irZ ^^^1 ~ TFT Pil ■ l^y) 



dz Hz ^'J'^'^' dz \Hz Hz 

In order to obtain an equation for the total B — L asymmetry Nb-Li we have flrst to 
write these two equations in the same flavour basis, for convenience the lepton flavour 
basis T-Tf*", and then subtract them. The rotation matrices are then given by 

for leptons and anti-leptons respectively, where notice that for matrices we indicate the 
heavy neutrino flavour index with an superscript in round brackets. Also notice that at 
tree level, corresponding to neglect CP violation, they would simply be conjugated of each 
other, i.e. 





In the charged lepton flavour basis one can flnally write the equation for the B — L 
asymmetry matrix as 

_ Dt(l) ^ ij d(1) D(l)t ^^^ij d(1) (on\ 



^For the sake of brevity in the fohowing we will say 'same (or different) flavour composition' and 'same 
(or different) flavour basis' for leptons and anti-leptons meaning that they are (are not) CP conjugated 
of each other. 
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whose trace gives the B — L asymmetry Nb-l- Notice that we are now identifying the 
a and a indexes since the lepton and anti-lepton flavor bases are CP conjugated of each 
other. We will find useful to define the following two matrices: 



V. 



(1) 



V. 



(1) 




b: 



(1) 



D(l)t 



(1) 



which, at tree level, are simply given by 



^o(i) _ pO(l)t 

' a/3 ~ ^ai 




Pi 
->0 nO* 



PlT 

-Clr Cl^i_ 

PlT 

— Clr C* , 

1 




(33) 
(34) 



no no* „0 



{hm)ii 




Using these results, we can now rewrite Eq. (32) as 



Di (iV^v, - iV^'J - Wi (pO(i)iV«-^; 



«/3 



dz 

The CP asymmetry matrix for the lightest RH neutrino A^^i is given by 



-(1) 

-a/3 



Vg} - Vi 
Ti + Vi 



(35) 



(36) 



(37) 



where the diagonal terms simply correspond to the flavoured CP asymmetries, Sal = Sia-, 
while the off-diagonal terms obey e^^l = (s^al)* and are not not necessarily real. This 



(i) 

expression can be generalised to the CP asymmetry matrix e)^j^, of any RH neutrino 
species Ni that in terms of the Yukawa couplings can be written as 



327r(/it/i). 



This expression slightly differs from that one in pTl [19] (simply the first term there is 
minus the imaginary part of the first term written here, so that the off-diagonal terms are 
real) while it coincides with the expression given in Notice that we have neglected 
terms 0{^p Nb-l) so that 




ROW ^0(l)t 



Nb-l 




ROW 



0(1) 



al3 



(39) 
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In the case we are discussing, when charged lepton interactions are neghg ible, and 
N^~^ commute. However, when charged lepton interactions become effective, and 
jsj-B-L commute in general and the term ('P'^(^) N^~^)ai3 is replaced by the anti- 

commutator til]: A^^-^)a/3 ^ {V°'^^\N^~^}^g/2. 

The diagonal components of the Eq. ( [32] ) can be explicitly written as 

" = e^^^D,{N^,-N^^\)-plW,Ns-L, (40) 



dz 
dz 



-(1) 



D 



B-L 



(41) 



Summing these two equations, one finally recovers the usual Eq. (jsj) for the total B — L 
asymmetry Nb-l = Tr[Ai'^o"^], which is washed out in the usual way at the production. 



On the other hand, from Eqs. (|40|) and (|4T|), one finds the relation 

1 



1 dN. 



Pi 



B-L 



dN 



B-L 



dz 



rP 



dz 



~^ 

Plr P^l 



Ni 



(42) 



which, together with Eq. ([2]), forms a system of equations that can be solved analytically. 
At low temperatures T <^ Tbi = Mi/zbi ^ Mi, the final values are then found to be 

Apir 



N. 



B-LJ 



P\rN'B-L 



N 



N 



B-L,{ 



„0 



N 



B-L 



+ 



2 



TVi ) 



(43) 
(44) 



This solution shows that the flavoured asymmetries contain terms that escape the washout 
at the production: these are the phantom terms [18]. If one only considers the unfiavoured 
regime, where charged lepton interactions can be neglected as we have done so far, the 
flavoured asymmetries are not themselves measured and the phantom terms cannot give 
any physical effect, in particular they cannot affect the baryon asymmetry. 

On the other hand, at T ~ 10^^ GeV, tauon lepton interactions start to be in equi- 
librium breaking the coherence of the lepton quantum states. These are described by 
additional terms in Eq. (32), which then generalises into pT] |3T| [25] 



dN^-^ 
dz 



-(1) 

-Q/3 



+ 1 



. Re(A,) 



Hz 

where we defined A'^, 





:) 











al3 



Im(A^) 

Hz 




(45) 




.N 



B-L 



al3 



al3 — ^^aP ^ ^"ap- 

self-energy are respectively given by 



The real and imaginary parts of the tau-lepton 
1331 



Re (A, 



64' 



-T 



and Im(A^) 
13 



8 X 10" 



(46) 



where is the tauon Yukawa couphng. The commutator structure in the third term on 



the RHS of Eq. (45) accounts for oscillations in flavor space driven by the real part of the 
self energy, and the double commutator accounts for damping of the off-diagonal terms 
driven by the imaginary part of the self energy. These effects can also be conveniently 
expressed in terms of the Pauli matrices cxj, i = 1,2, 3, so that the evolution equation 
reduces to 

(IAtB-L ^ 

= 4>i(iVM-iV^^J-^H^i{^°^^\iV^-^U (47) 

In order to close the system of equations, we also need a Boltzmann equation for the A^£+£, 
which is given by 

^ = -^jf^i-2UN^i' - S, (Nif - 2 , (49) 

where Sg = Tg/{Hz) accounts for gauge interactions. As shown in [53], this term has the 
effect of damping the flavor oscillations. This can be understood by noticing that gauge 
interactions force N^^^ = 2N^^6a/s, which in turn renders the oscillatory term Eq. (45) 
negligible. 

Let us indicate with ^ 10^^ GeV that value of the temperature below which one can 
approximate, with the desired precision, the lepton quantum states as a fully incoherent 
mixture of |r) and \t^) quantum states corresponding to a complete damping of the off- 
diagonal terms in the lepton density matrix (analogosly for anti-leptons) . This means 
that for T < the r and lepton asymmetries, given by the diagonal entries of 
N^p^, are fully measured by the thermal bath and reprocessed by sphaleron processes 
conserving the and A^x asymmetries, so that at T* one has A^J* = N!^~^[zb) and 
iVT* = N^-^^{zb). However, notice that since, for Mi > 10^^ QgV, the total B - L 
asymmetry has been already produced and got frozen in the unfiavoured regime, this 
fully flavoured regime stage does not affect the final total B — L asymmetry. Even with 
-^TVi 7^ 0) phantom terms do not contribute to the final asymmetry because they cancel 
with each other. In other words, within the A^i-dominated scenario, phantom terms have 
no effect on the final asymmetry]^ Phantom terms are therefore a consequence of charged 



^On the other hand remember that if one considers 10^ GeV ^ Mi <^ 10^^ GeV, the two fully 
flavoured regime holds during the period of leptogenesis and the density matrix equations reduce to the 



set of classical Boltzmann equations Eqs. (16). The terms in the flavoured asymmetries coming from 



CP violating terms due to a different flavour composition are still present but they are not phantom, 
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lepton flavour effects but they do not have any consequence on the final asymmetry in 
the A'^i-dominated scenario in the absence of heavy neutrino flavour effects. 

Our results show explicitly the presence of the phantom terms. They extend previous 
results [ini EH] where the lepton and anti-lepton density matrices were, implicitly or 
explicitly, assumed to be diagonalizable in bases that are CP conjugated of each other and 
this precludes the derivation of the phantom terms. However, as far as the A'^i-dominated 
scenario is considered, phantom terms can be safely neglected in the calculation of the 
final asymmetry, and therefore there is no contradiction between our results and those of 

da Eg. 

On the other hand, since we are interested in accounting for heavy neutrino flavour 
effects, we cannot neglect phantom terms. Indeed, as we will discuss, they can in this case 
affect the total final asymmetry giving rise to the possibility of phantom leptogenesis [18] . 

It is also quite interesting to notice that, even in the A'^i-dominated scenario, the 
phantom terms can produce lepton flavour asymmetries much larger than the observed 
baryon asymmetry, though depending on the initial RH neutrino abundances. These 
large cancelling flavour asymmetries can be as large as 10~^ — 10~^ for initial thermal 
Ni abundance and in this case they could have potential effects at lower temperatures 
T < 100 MeV when the thermal bath is flavour sensitive. They can very unlikely have 
direct effects on primordial nuclear abundances since an impact on BBN requires much 
larger asymmetries (> 0.01 [M])- However, asymmetries as small as ~ 10^^ could for 
example be relevant for active-sterile neutrino oscillations [35]. 



3 A simplified case with two charged lepton fiavours 

In this section we account for heavy neutrino flavour effects considering a simplifled two 
charged lepton flavour case. This will greatly simplify the notation making the new 
results more easily readable. It will be then quite straightforward in the next Section to 
generalise all the equations to a realistic three lepton flavour case. For deflniteness we 
consider masses Mi ^ 10^ GeV, when only tauon lepton interactions have to be taken 
into account. 

We also assume that the heaviest RH neutrinos N^, do not contribute to the flnal 
asymmetry. This is in any case a valid assumption if M3 ^ Tjih ^ M2, since in this way 



since they are measured directly at production and undergo washout. Therefore, if there is a flavour- 



asymmetric production, they contribute to the final asymmetry, yielding the second term in the Eq. (171, 
and can even dominate as we discussed. 
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Figure 3: Flavour configuration of tlie two lieavy neutrino lepton flavours, and ^2■, lead- 
ing to the simplified two charged lepton neutrino flavour case considered in this section. 



the A^a's would not thermalise. As for lepton flavours, we will extend the results to the 
three heavy neutrino flavour case in the next Section. 

Notice that with these assumptions, the two charged lepton flavour case can be re- 
garded as a special case where the two heavy neutrino lepton flavours, £i and £2, lie on 
the same plane orthogonal to the e — fi plane and therefore = = t"*" (see Fig. 3). 
Correspondingly the two anti- lepton flavours, ii and £2, also lie on the same plane or- 



thogonal to the e — /i plane and therefore T2 



with r-*- that is now assumed to 



be CP conjugated of r"*". In this way, in the whole following discussion, we will have only 
two charged lepton flavours, r and r-*-. 



The density matrix equation Eq. (|45|) for N^n ^ , valid for A^i leptogenesis, gets then 



a/3 



generalised into 



dz 



-(1) 



Ni 



+ i Re (A, 



a/3 



(50) 





:) 











- Im(A^) 



a/3 





.A^ 



B-L 



We now discuss the three asymptotic cases, the first one for M2 ^ 10^^ GeV ^ Mi, the 
second for M2, Mi > 10^^ GeV and the third for Mi, M2 < 10^^ GeV, where Boltzmann 
equations are recovered. In this way we will derive, within a density matrix formalism. 
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results that were already obtained within an instantaneous collapse of the quantum state 
formalism: the first is phantom leptogenesis [TH] , the second is the heavy neutrino flavour 
projection [6l 120] . 



3.1 Case M2 ^ 10^^ GeV ^ Mi: three stages phantom leptogene- 
sis 

Let us consider the asymmetry production from the N2^s at T ~ M2. This is basically de- 
scribed by the same equations that we wrote in the previous section for the A'^i-dominated 
scenario where now simply all quantities have to be relabelled in a way that 1 — )■ 2. Since 



charged lepton interactions are negligible, we can use the Eq. (36) for the calculation of 
^aB^i ^^^t '^o^ simply becomes 



-(2) 



D2 {N, 



N2 



(51) 



dz 

with an obvious re-definition of all quantities that now refer to A^2- At T ~ Tb2 = ^ij zbi 



the r and r asymmetries are described by the Eqs. (43) with 1 — )■ 2 



Tbi) 



P2 



j\^T~Tb2 



2 



^^N2 ' 



'B-L 



+ 



N 



N2 



(52) 
(53) 



Again, at temperatures below T* <^ 10^^ GeV, the iV^"^ off-diagonal terms have been 
fully damped by the tauon charged interactions, so that the A^a^ and A'a^x asymmetries, 
corresponding to the diagonal terms, can be treated as measured quantities. 

At T ~ Tb2, the phantom terms in the Eqs. (52) cancel with each other and they do 



not contribute to the total asymmetry. Therefore, so far, the description of the asymmetry 
evolution is completely analogous to that one in the A^i-dominated scenario. 

However, there is still a third stage to be taken into account: the lightest RH neutrino 
washout. For T ~ Mi, the tauon and the r-*- asymmetries are washed out by the lightest 
RH neutrino inverse processes. At T ~ Tbi = Mi/zbi, they get frozen to their final 
values 



P2t ^^B-L 



'2r 



N2 



^0 atT~Tb2 _j_ ^P2r ^in 



B-L 



N2 



37r Ty- 



(54) 
(55) 
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so that the final total asymmetry N^b-l — ^At + ^At^' fiavour a = t (r-*-) one 

has ^ 1, while for the other fiavour /3 = (r) one has Ki^ ^ 1, the final asymmetry 
will be dominated by the a asymmetry, 

NL,^pLnIzI^^-^N-. (56) 

The interesting thing is that the phantom term is not affected by the washout at the 
production, in a way that, if this is very strong, the phantom term can dominate the 
final asymmetry if A^^^ is non- vanishing and sufficiently large (phantom leptogenesis) . 
This result has different implications. As noticed [HI [IT], phantom terms introduce a 
further source of dependence of the final asymmetry on the initial conditions. Therefore, 
phantom leptogenesis scenarios necessary rely on a description of the initial conditions. 
On the other hand, as we already noticed in the case of the A'^i-dominated scenario, large 
phantom terms could also have some relevance at lower temperatures when the thermal 
bath is fiavour sensitive. 

Phantom leptogenesis has been first discussed within an instantaneous quantum state 
collapse description [TH], here we have re-derived it within a density matrix formalism. 
Notice that there are three well separated stages: N2 asymmetry production at T ~ Tb2, 
decoherence at T ~ and asymmetric A^i washout at T ~ Mi. 

Notice also that phantom leptogenesis has some analogies with the scenario of A^^i- 
leptogenesis with £1 = [10] that we discussed in the previous section. In both cases the 
final asymmetry originate from the CP violating terms oc Api^ due to a different flavour 
composition of leptons and anti-leptons. In both cases a non-vanishing flnal asymmetry 
relies on a asymmetric washout acting on the two flavour asymmetries. There are however 
important differences. In the case of A'^i leptogenesis with £1 = one has that production, 
decoherence and washout occur simultaneously, while in the case of phantom leptogenesis 
they occur at different stages and between the production and the Ni washout stage 
the phantom terms they cancel in the flnal asymmetry. Another important difference is 
that in the case of phantom leptogenesis one has not to assume the special assumption 
El or €2 = {B — L conservation): it is enough that the washout at the production 
is sufficiently strong that phantom terms can naturally dominate. On the other hand 
phantom leptogenesis rely on non-vanishing and sufficiently large initial N2 abundances. 

As we are going to show, phantom leptogenesis is even more general and it does 
not necessarily require that the N2 production and the A''i washout stages occur in two 
different fully flavoured regimes. 
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3.2 Case M2 > 3 Mi ^ 10 GeV: heavy neutrino flavour projec- 
tion and two stages phantom leptogenesis 

Let us now consider the case when both heavy neutrino masses M2, Mi ^ 10^^ GeV and 
charged lepton interactions do not affect the final asymmetry. This can be called the heavy 
flavoured scenario jl^ since the only lepton flavours that affect the final asymmetry are 
those produced from the heavy RH neutrinos. 



3.2.1 Projection effect (in isolation) 

For illustrative purposes, we want first to describe just the A^^i washout of the asymmetry 
produced by the N2 decays, without any additional effect. Therefore, we first neglect the 
different flavour compositions of leptons and anti-leptons assuming that Apia = Ap2a = 0. 
With such a simplifying assumption, the lepton quantum states are given by 

|1) = Cir \r) + Ci.x |r^) and |1) = ^ |f ) + C^V |r^) , (57) 

\2)=C2r\T)+C2r^\r^) aud \2)=Cl\f)+C'^-^\f^). (58) 

Assuming the hierarchical limit, M2 > 3 Mi [36], there are two well distinguished different 
stages. In a first stage at T ~ M2, an asymmetry is produced from N2 decays. The lepton 
density matrix is then given by pfj = diag(l, 0) in the basis £2 — ^2- Analogously the anti- 
lepton density matrix is given by p|j = diag(l, 0) in the basis £2 — ^2 that for the moment 
we are assuming to be CP conjugated of £2 — £2- in the previous subsection, the 



asymmetry production from N2 decays is again described by the Eq. (51) with vanishing 
phantom terms so that we simply have 

Arf,-^(r ^ Tb2) ^ pI , nX'L {T ^ Ts2) ^ vl^ AT^I?^ , (59) 

where N'^Zif^ ~ e2K,{K2). We have now to consider the Ni washout stage at T ~ Mi. 
Since at the moment we are just interested in describing the A^^i washout, we also neglect 
the A^^i asymmetry production assuming a vanishing e^J^^. Moreover let us first further 
assume |1) = |r) and correspondingly |1) = |f). 



In this way, at T ~ Mi, the Eq. (50) for the asymmetry evolution in the charged 



lepton flavour basis simply reduces to (a, (3 = t,t^ 

and, at the end of the A'^i-washout at T 
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Finally, at T ~ 10^^ GeV, the charged lepton interactions damp the off-diagonal terms 
measuring the tauon and the 'non-tauon' (i.e. the r-*-) asymmetries. 

This result can be easily generalised. Let us, first of all, allow an arbitrary |1) flavour 
composition but continuing, for the time being, to neglect the Ni asymmetry production, 
at T ~ Tei. The Eq. (60) has now to be written in the basis ii — ii, 



dz M iN^'r^ 



I ' r I (^i,Ji = l,l")- (62) 



The solution is again quite trivial in this basis: the 11 term is washed out, 

iVff^(T ^ Tbi) = e-^ Nt\T ^ Tb2) , (63) 

while all the other terms are unwashed. The asymmetry matrix at T ~ Tb2, in the ii — i^ 
basis, can now be calculated in terms of the rotation matrices (cf. Eq. (31)) as 

A'--(T Ts,) = Nlzl"' C C ( J ° ) ■ (64) 

In a more compact way, considering that N^~^{T ^ Tb2) = ^B-r^^|2)(2|, this can be 
more conveniently written as 

ATB-L(rj.^rp N ArT-TB2 ( ^'12 (1|2)(2|1^) \ 

iV„,, iT^Ts,) = N,_, ((j,|2)(2n^ j , (65) 

where [IT] 

The final asymmetry can then be calculated as 

N'b-l = Tr[iV^^^(T ^ T^O] = e"^ P12 + (1 - P12) . (67) 

The asymmetry can be also rotated in the charged lepton flavour basis, 

N^^^{T ^ Tb,) = Rlf Nf^i^^T ^ Tb,) R^J . (68) 

At T ^ 10^^ GeV the charged lepton interactions just damp the off-diagonal terms without 
affecting the total asymmetry given by the trace and for this reason we could directly write 
the Eq. 

This result fully conflrms what one expects within an instantaneous quantum state 
collapse description. It is not only conflrmed that just the orthogonal component (in our 
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simplified case the r-*- component) undergoes the A^i washout while the parallel component 
(the r component) completely escapes it [21120], but also that the washout of the parallel 
component is exactly described by the factor exp[— (Svri^i/S)], independently of the value 
of Ki [17]. Notice that in an intermediate regime /^i ~ 1, the quantum states at T ~ Tei 
are left in a sort of partially incoherent mixture, with some residual flavour oscillations 
that however do not affect the total asymmetry. 

Notice that this result also applies to a possible pre-existing asymmetry produced by 
some other external mechanism [201 E]- Therefore, the conclusions of [17], employing 
this result in various situations, are also confirmed. 

One can then easily further generalise this result accounting also for a possible A^i 
asymmetry generation, simply obtaining for the final asymmetry 



£i«(Ki)+ (e-^^'^Pi2 + l-pi2) 62 ^{K2). (69) 



3.2.2 Projection effect in combination with phantom leptogenesis 

We still miss the (usual) last step. We have so far assumed that the flavour compositions of 
the £2 and £2 quantum states are the same. We want now to show that when this additional 
flavoured CP violation contribution is taken into account, again phantom terms appear 



in the expression Eq. (69) 



Therefore, the quantum states |2) and |2) have now to be written more generally as 

|2) = (1|2)|1) + (1^|2)|1^) and |2) = (1|2) |i) + (1^|2) |1^) . (70) 

Notice that this time the role played by the charged lepton flavour basis in the previous 
subsection, is replaced by the heavy neutrino lepton basis £i-£i ■ Therefore, writing the 



Eq. (51) in this basis, we obtain 



-^-e.,,,D2{N^2-N^^)-W2Ns-L ^ ^iX|2)0^2|l)0 1 - rf, j ' ^^'^ 

where the superscript "0" indicates, as usual, the tree level quantities that can be approxi- 
mately employed in the calculation of the washout term. In this way we obtain expressions 



for the heavy neutrino flavoured asymmetries, that are analogous to the Eqs. (52) for the 
charged lepton flavoured asymmetries. 



N^,-\Tc^Tb2) ^ pl2e2K{K2)-^N'^^, (72) 
iVf,-^(T^T^2) ^ {l-p%)e2K{K2) + ^N^^^. (73) 
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Finally, taking into account the lightest RH neutrino washout, we obtain for the final 
asymmetry 

N'^_^=e,K{K,) + 

Therefore, the phantom terms give an additional contribution to both components and 
in particular to the orthogonal component. As usual this contribution is (A^2) unwashed 
at the production and proportional to the N2 initial abundance. \{ Ki ^ 1, both the 
parallel and the orthogonal components are unwashed and the phantom terms cancel with 
each other. On the other hand, in the opposite case, for i^i ^ 1, the parallel component 
is completely washed out so that only the orthogonal one survives together with the 
additional (A^^i) unwashed phantom term contribution. 

This result shows that phantom leptogenesis goes even beyond the case where the two 
RH neutrinos masses fall into two different fiavour regimes [18]. 

Finally, it should be clear that an account of the different flavour compositions of the 
^l and C-i quantum states, leads to additional phantom terms. These, however, cancel 
with each other and do not contribute to the flnal asymmetry, as already discussed in 
section 2. 



Sir 
8 



37r TV' 



12 



(74) 



3.3 Case lO^^ Q^y > Mi, M2 

When Ml, M2 ^ 10^^ GeV both RH neutrinos produce their asymmetry in the two-flavour 
regime. The production from the heavier RH neutrino is given by the usual result 

N^-\T ~ Tb2) = e2r '^{K2r) , N^^T :^ Tb2) = e2r^ ^{K2r^) . (75) 

In the strong washout regime for both flavours, K2t^iK2t ^ 1, the sum, the total asym- 
metry, can be approximated by the Eq. (17) rewritten for the heavier RH neutrino. When 
the temperature drops down to T ~ Tgi, the washout from the lighter RH neutrino starts 
to act. Similarly to the previous cases, this washout factorizes from the general expression 
and can be expressed as a simple exponential pre-factor. We have 

N^-\T^Tb2) = £2./t(i^2.)e-^^-, (76) 
NXiiT^TB2) = e2r^<K2r.)e-'f^^^-. (77) 

The production of the asymmetry from the A^^i decays is then added to what is left from 
the production, so that we flnally obtain 

N^-\T^Tbi) = e2rK{K2r)e-'^'''^ +eirK{Kir), (78) 
NX^{T^Tbi) = e2r^ K{K2r^)e~'^''^r^ +6,,^ k{K,,^). (79) 
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Figure 4: A generic three heavy neutrino lepton flavours configuration. 



It should be noted that there are no phantom terms in this formula because of the as- 
sumption made at the beginning of the Section that = ^i" = t-*-. In this case, we 
have an effective two-fiavour problem and there is no hidden asymmetry left unwashed 
at production and, therefore, no phantom terms. If we relax this assumption and have 
7^ Tf*", we have to work in a full three-flavour basis, and, as we will see, phantom terms 
will appear again in the flnal asymmetry. We will study such a case in the next Section. 

4 General case with three charged lepton flavours 
and three heavy neutrino flavours 

If we consider the general realistic case where the full three lepton flavour space is con- 
sidered, the density matrix equation has to be written in term of 3 x 3 matrices. In 
general the three heavy neutrino flavours have no particular flavour orientations in the 
three charged lepton flavour space (see Fig. 4). If we also consider generic three RH 



neutrinos mass patterns with masses Mj ^ 10^ GeV, the density matrix equations (50) 
further generalise into (a, (3 = t, fi, e) 
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\W,{V'^^\N^-^}^^ (80) 
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If one of the three masses is lower than ~ 10® GeV, electron flavour interactions terms 
have to be included as well, though they have no real impact, within this framework, 
on the final asymmetry. This is because the electron asymmetry is in any case already 
measured as a 'neither-muon-nor-tauon' asymmetry. 

This master equation can now be used to calculate the final asymmetry not only for 
all the 10 mass patterns shown in Fig. 1, but also when the Mj's fall in one of the flavour 
transition regimes. 

Notice that, though in this paper we are only considering hierarchical RH neutrino 
mass patterns, this equation can also be used to calculate the asymmetry beyond the 
hierarchical limit |36j and even in the resonant case [37] . In this latter case, however, 
many different effects can become important and should be included |38j . 

Solutions of this set of equations are particularly difficult when at least two of the 
five kinds of interactions are simultaneously effective, something that goes beyond our 
objectives. Here, as an example with three flavours, we want to show a particularly 
interesting asymptotic limit that could not be described within the simplified two-flavour 
case discussed in the previous section: the two RH neutrino model [15]. We will show 
that, even in this case, phantom terms have in general to be taken into account. 



dz 



4>i (Nn.-KV- 



+ 4>2(iV/N 



A^' 



oq - 



+ 4>3 (iV/v3 - iV^l) 



4.1 Boltzmann equations for the two RH neutrino model 

We consider a two RH neutrino model [39jJ corresponding to a situation where M3 is 
sufficiently large (M3 ^ 10^^ GeV) to decouple in the seesaw formula for the calculation 
of the neutrino masses |40j . In order to reproduce the observed baryon asymmetry one has 



to impose Mi > 10 GeV so that the muon interactions can be neglected in the Eq. (80) 
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On the other hand, in order to have Mi and M2 as low as possible, it is interesting to 
consider the case 10^ GeV ^ Mi <^ M2 ^ 10^^ GeV in a way to obtain a RH neutrino 
mass spectrum corresponding to the third panel (starting from upper left) in Fig. 1. 
This model has been recently revisited in [15]. We want here to re-derive, starting from 



the density matrix equation Eq. (80), the Boltzmann kinetic equations and the consequent 
formula for the final asymmetry that in [T3] has been used to calculate the value of Mi 
necessary to reproduce the observed baryon asymmetry |^ 

Thanks to the hierarchical limit, we can again introduce different simplifications. First 
of all we can impose the complete damping of the ra and ar {a 7^ r) off-diagonal terms 
in the asymmetry matrix. 

Second, we can consider the N2 production at T ~ Tb2- With these assumptions, only 



the A^2-terms can be considered in the Eq. (80) and the asymmetry matrix can be treated 



as a 2 X 2 matrix in r-r^ flavour space. In this way the density matrix equation reduce 
to a set of two Boltzmann equations in an effective full two flavour regime. 



dN. 



B-L 

TT 



dz 

dN^'r'^ 



dz 



ef} D2 (iV^, - iV^^) - pI W2 Nf-'^ , (81) 
't^^ D2 (Nn. -KV~ <x W2 <4 . (82) 



As usual, assuming first for simplicity that the |r^) and the |r^) quantum states have 
the same flavour compositions, one finds 



N^-L(^T^TB2)=e2rK{K2r) and NX'^{T^TB2) = e2.^K{K2r^). (83) 



These values of the asymmetries at the end of the N2 production stage have to be used as 
initial values in the set of equations describing the evolution of the asymmetries during 
the Ni production. 



dN. 



B-L 

TT 



dz 
dN^r"! 



dz 



s^t} Di {N^, - N^^\ ) - pI Wi N^^-' , (84) 
e^X^ Di {Nr,, - ) - p° Wi N^i ■ (85) 



The component of the asymmetry at the end of the N2 production has to be decom- 
posed into a parallel component and into a orthogonal component that we indicate 
with the symbol r^^. In this way one finds that the final asymmetry is the sum of three 
flavour components (see Fig. 5), 



^It has been shown in [T5] that even the N2 production depends just on Mi and not on M2, provided 
that this is much smaher than 10^^ GeV. 
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Figure 5: Relevant lepton flavours in the two RH neutrino model. 



N'j,_^ = iV«-^(T ^ Tbi) + iVf,;i(T ^ Tbi) + iVf,-^ (T ^ Tbi) , (86) 

where 

iVf^-^(T ~ Tbi) = £l./^(irl.)+£2r/^(i^2r)e-^^l^ (87) 

N^f^^AT^TBi) = e,,.K{K,,.)+pl^^^^.e,,.K{K,,.)e-'^''^^^\ (88) 

iV5-;^x (T ^ Tsi) = (l-p\^^e,,.K{K,,.). (89) 

This expression coincides with the result found in [15] and is strictly valid for vanishing 
initial N2 abundance. 

More generally, if N^^ 7^ 0, one has also to include the phantom terms allowing 
for different flavour compositions between the |r^) and the |f^) quantum states. The 
procedure is essentially the same as in Section 3.2, with the only difference that now the 
phantom terms will appear only in the and t^^ components but not in the measured 
tauon component. We can therefore directly write the final result, 

N'b_l = Sir K{K,r) + £2r f^iK2r) e^'-i (90) 

(An ± ± \ _32L A- 
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where each of the three hnes corresponds respectively to the r, rf^ and t:^^ components. 
This last example shows, once more, how phantom terms are present all times that the 
production occurs either in a two or in a three flavour regime though only those generated 
by the heavier RH neutrinos can be afterwards asymmetrically washed out by the lighter 
RH neutrinos and contribute to the final asymmetry without cancelling with each other. 

5 Final discussion 

Within a Boltzmann classical kinetic formalism one has to distinguish the ten different 
RH neutrino mass patterns shown in Fig. 1. These are obtained in the limits where the 
masses Mi are hierarchical and do not fall in the transition regimes. We have extended 
the density matrix formalism for the calculation of the matter-anti matter asymmetry in 
leptogenesis including heavy neutrino flavours. In this way we obtained a unique density 
matrix equation describing the asymmetry for any choice of the RH neutrino masses, even 
beyond the hierarchical limit. 

Within this more general description, the ten hierarchical RH neutrino mass patterns 
of Fig. 1 correspond to those cases where the (five) different interactions are only singularly 
effective within a given range of temperatures. In this way the evolution of the asymmetry 
can be described in well separated stages where the density matrix equations greatly 
simplify reducing to multiple sets of Boltzmann equations, one for each stage. In this 
cases we recovered and extended results that were already derived within a description 
based on an instantaneous collapse of lepton quantum states. 

The flavour projection effect, where the orthogonal component of a previously pro- 
duced asymmetry escapes the RH neutrino washout, is fully confirmed. We have also 
shown that the washout of the parallel component is exactly described by the usual ex- 
ponential washout factor independently of the washout regime. 

Phantom terms emerge as quite a generic feature of flavoured leptogenesis and have 
to be taken into account for non-vanishing initial RH neutrino abundances. We have 
indeed shown how their presence goes beyond the A''2-dominated scenario where they 
were originally discussed [18]. They can contribute to the flnal asymmetry even if the 
production from an heavier RH neutrino species and the washout from a lighter RH 
neutrino species occur in the same flavour regime. We have also noticed that, even though 
in the iVi-dominated scenario they exactly cancel with each other and do not contribute 
to the flnal asymmetry, they can give rise to flavoured asymmetries much larger than 
the observed baryon asymmetry with potential interesting applications, for example in 
active-sterile neutrino oscillations. 
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Even though we have exphcitly calculated the final asymmetry only in one of the ten 
asymptotic limits RH neutrino mass patterns shown in Fig. 1, in the case of the two 
RH neutrino model, the procedure can be easily extended to all others neutrino mass 
patterns. For example one can easily show the expression for the final asymmetry in the 
iV2-dominated scenario, when Mi < 10^ Gev [T3l fTHIiH [T8] . 

It would be desirable in future to calculate the asymmetry beyond these ten asymp- 
totic limits, solving the full density matrix equation. In this way the calculation of the 
matter-anti matter asymmetry would be extended to a generic RH neutrino mass pattern, 
including the cases where the RH neutrino masses fall in the transition regimes between 
two fully fiavoured regimes and quantum decoherence from charged lepton interactions 
acts simultaneously with asymmetry generation and wash-out. This would make possi- 
ble to interpolate between the asymptotic limits, finding the exact conditions on the RH 
neutrino masses for the validity of the solutions that we have discussed here. 
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